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Summary. — 

The electron and positron cosmic rays puzzle has triggered a revolution in the 
field of astroparticle physics. Many hypotheses have been proposed to explain 
the unexpected rise of the positron fraction, observed by HEAT and PAMELA 
experiments, for energies larger than few GeVs. In this work, we study sources 
of positron cosmic rays related to annihilation of dark matter and secondary 
production. In both cases, we consider the impact of uncertainties related to 
dark matter physics, nuclear physics and propagation of cosmic rays, finding that 
the largest uncertainties come from propagation. We find that some key-features 
present in the positron signal from dark matter annihilation are preserved even 
though the uncertainties. In addition, we study the stability of the positron fraction 
under small variations of the electron fiux, which is usually considered as known, 
we found that considering just the observational uncertainties in the electron flux is 
enough to changed dramatically the positron fraction in the energy range were the 
excess is observed. 
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1. Introduction 

During the last decades, the astrophysical and cosmological evidences of dark mat- 
ter and dark energy have created a revolution in the field of fundamental physics. Some 
models of new physics, which aim to extend the Standard Model of particles, predict par- 
ticles vifith the right properties to act as dark matter candidates. In a similar way, cosmic 
ray physics is continuously stimulated with new observations showing an environment 
much more active than it was expected. Furthermore, cosmic rays correspond to genuine 
samples of the matter composition of the Milky Way and give essential information about 
the Interstellar medium and the Solar System environment. 
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A fraction of comic rays is composed by antimatter. This component may provide 
crucial clues on the non-standard sources of cosmic rays because it is less abundant than 
the matter cosmic-rays component. This analysis may open a window for dark matter 
searches. 

This paper is based on results of the author's Ph.D. Thesis [1] and some other projects 
derived from it. In Section 2, we study the standard mechanisms of production of elec- 
trons and positrons. In the same spirit, in Section 3, the propagation of cosmic rays is 
reviewed and discussed. Then, we study the production of positrons in scenarios of dark 
matter annihilation (Section 4) and secondary production (Section 5). 



2. Production of positrons and electrons 

In a first approach, the Standard Model of particles physics describes well enough the 
physics behind the production of cosmic rays. However, let us clarify that, the produc- 
tion in the astrophysical context arises both from particles physics and from the galactic 
environment. 

Among the sources of cosmic rays, supernovae are the main ones and a big fraction 
of them comes from the matter expelled at the time of the explosion. This mechanism 
is generally the dominant one but not for all the species and energy scales. In fact, 
interactions between cosmic rays and the interstellar gas sizeably contributes, adding an 
extra amount of cosmic rays known as secondaries. 

Electrons and positrons are also produced by other mechanisms, like pulsar emis- 
sion [2], and decay /annihilation of dark matter [3]. 



2'1. Production in proton-proton collisions. - The nuclear scattering is the 
fundamental process behind secondary electrons and positrons. The secodary production 
takes place in the galactic disk when nuclei cosmic rays scatter off the interstellar medium. 
In particular, it depends on the composition of the interstellar medium - composed by 
hydrogen (0.9 part. /cm^) and helium gas (0.1 part. /cm"^) [4, 5, 6, 7] - and on the fluxes 
of protons and alpha particles [8]. 

2'1.1. Production from mesons decay. In proton-proton collisions, charged pions 
and kaons are produced in big amounts. Their decay channels are responsible of the big 
amount of electrons and positrons. Of course, this is valid for collisions at low energy 
(piab < 20GeV) when meson production is dominant. 

The electron and positron energy spectra are calculated from the inclusive cross sec- 
tion of pions and kaons and from the electron and positron spectra associated to their 
decays: 

dSe dse dee ' 
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where a particular contribution is generically calculated as: 

2) — (£e,£p = / dex ' £x,£p X /e.x(£e,£x) , 

d£e J dex 

where X denotes pions or kaons as intermediate particles and fe,x the meson decay dis- 
tribution into electrons or positrons. The decay distributions is analytically calculated 
including muon polarization effects [1]. Pion decays are the simplest ones because they 
have just one dominant decay mode (/x^z^^) instead of kaon decays which are complex to 
treat (details in appendices C and D of [1]). 

The inclusive cross sections of pions and kaons from proton-proton collisions are very 
difficult to calculate a •priori. In the literature, two well known parametrization for the 
inclusive cross section of pions and kaons are available: 

• Badhwar et al. parametrization [9, 10]. 

• Tan-Ng parametrization [11, 12]. 

Both parametrization reproduce good enough the observations in colliders. In any case, 
the small differences among them gives an estimation of uncertainties coming from nu- 
clear and particle physics. 

2'1.2. Kamae et al. PARAMETRIZATION. A parametrization of production of positrons, 
electrons and other particles is proposed by Kamae et al. [13, 14, 15]. It aims to provide 
an easy way to compute and estimate cosmic ray fluxes that comes from interstellar 
medium interactions with nuclei cosmic rays. 

New processes are included like contributions from A(1238) and many hadron reso- 
nances around 1600 [MeV/c^] that make it accurate in the very low proton energy range 
(piab < 8 GeV) . Other included processes come from diffractive dissociation which con- 
tributes in the intermediate energy range {p\ah > 20 GeV) . 

In addition, this parametrization works for very high energy scales being based on 
simulators like PYTHIA [16] to generate the spectra at higher energies. Also, its big 
advantage is to save CPU time by avoiding to calculate the convolution among mesons 
production and meson decay into electrons and positrons. 

2'1.3. Uncertainties from nuclear physics. Previously we described different 
approaches to calculate the electron and positron spectra in proton-proton collisions. 
Each parametrization is based on physical assumptions and encompasses available ex- 
perimental data: therefore the small differences among them can be used to quantify the 
nuclear physics uncertainties. 

In Figure 1, we show the comparisons among the three parametrizations for different 
proton kinetic energies. We note that these parametrization are closely similar in behav- 
ior. However, there are variations up to 80% at proton energies of 20 GeV, as in case of 
Kamae et al. versus Tan and Ng parametrization. Another feature is that the Kamae's 
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Figure 1. ~ Positron and electron inclusive cross section versus energy for proton energies of 
2, 10 and 35 GeV. The different panels show the case of Badhwar et al. [9], Tan and Ng [12], 
and Kamae et al. [14] parametrizations. The last panel shows that the relative difference with 
respect to Tan and Ng solution is on average around 15% and 25% for Badhwar et al. and 
Kamae et al. cases respectively. 



paramctrization estimates a smaller electron cross section with respect to the other two 
parametrizations. 

Due to the low statistics at very low energy, Badhwar's and Tan's paramctrization 
tend to produce non-physical distributions for proton kinetic energies below 6 GeV. Nev- 
ertheless, the total inclusive cross section - the integrated version of those - are still in 
agreement with the available experimental data. To fix this undesirables feature, both 
paramctrization are patched by doing a smooth transition from 3 GeV until 7 GeV with 
the Steeker's model [17, 18]. Let's clarify that Kamae's paramctrization also includes 
that feature, but considering more resonances. 

Moreover, for proton energies above 100 GeV, Badhwar's paramctrization becomes un- 
stable specially the electron cross section case. 
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Table I. - Intermediate states for positron and electron multiplicity distributions. 



2'2. Production in annihilation processes. ~ Electron and positrons can also 
be the result of annihilations. A first case to come in mind is the matter-antimatter 
annihilation like electron positron collisions at LEP. 

The annihilation of dark matter also enters in this category, providing a new type of 
source of electrons and positrons that it would coexist with standard cosmic rays sources. 

Independently of how these particles annihilate, electrons and positrons would come 
directly from the annihilation event (direct production case) or from the annihilation's 
sub-products, like decay of gauge or higgs bosons or hadronization/decay of quarks. 

Wc work in a general approach, in which we generate electron and positron multiplic- 
ity distributions which are independent of dark matter physics and can be used for any 
purpose. 



2'2.1. Multiplicity distribution. In a generic annihilation case, the multiplicity 
distribution of electrons and positron is: 




where we decompose the annihilation in intermediate states i, related directly to the 
annihilation mechanism via branching ratios, for instance, the dark matter annihilation 
case: 

(4) BR(xx-.)^^i^^. 

Ctotal 

This presents a general decomposition based on bricks which are the multiplicity distri- 
bution for many states i. 



2 "2. 2. Calculation of Multiplicity distributions. We consider the most gen- 
eral set which is composed by pairs of Standard Model particles. Each pair preserve 
quantum numbers of the annihilation, i.e. an electrically neutral and colorless final state. 

Depending on the case, we calculate analytical distributions for simple cases (elec- 
trons, positrons, and muons), and for complex ones, like the ones which involves hadroniza- 
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Figure 2. - Multiplicity distribution for positrons and electrons versus energy. Each panel shows 
a diflterent intermediate state - W^W~ , t^t^ , bb, and tt - with dark matter masses (m^) of 
100 GeV (or 200 GeV) and 1 TeV. The different shapes on electron and positron distribution 
come from the effect of production of polarized muon. 



tion processes, we used a modified version of PYTHIA [16] to generate the distribution, 
which includes the efFect of polarization of muons. 

The basic set is based on Standard Model particles Table 1. Special cases like the 
Standard Model higgs and models like Two Higgs Doublet Models are easily composed 
using the basic set (details in [1]). 

In Figure 2, we present multiplicity distributions for cases inspired by annihilation 
of dark matter. We observe that the shape of each distribution depends directly on the 
intermediate state: states involving quarks produce a softer electron and positron spectra 
than leptonic cases. In addition, the effect related to the polarization of muon in meson 
decays also has an impact, producing more energetic electrons than positrons. 
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Figure 3. - Propagation Zone geometry for the Milky Way. A cylinder of radius R (= 20 kpc) 
with thickness of 2 L delimits the region where cosmic rays propagate. A small cylinder with 
same radius but with thickness 2 {— 200 pc) models the Galactic plane. The Solar system is 
placed at the Galactic plane with distance tq = 8.5 kpc from the Galactic center. 

3. Propagation of positron and electron cosmic rays 

The cosmic rays' journey from the source until the Earth is a complex problem. Since 
cosmic rays start to travel, they are affected by many processes and their intensity de- 
pends on the cosmic ray energy scale. For instance, in the GeV-range, their are strongly 
affected by turbulent magnetic fields that induce spatial diffusion, similar to the behavior 
of a single molecule in a gas. 

For the GeV-scale cosmic rays, the Galactic medium has a very important role on the 
propagation. Apart from the turbulent magnetic fields, there are also interactions with 
diffuse radiation fields (UV, IR, CMB). Their continuous interaction changes the cosmic 
rays energy, electrons and positrons are more affected than other species. 
The combined effect of diffusion, energy losses, distance, among others, cause that the 
observed cosmic rays spectra is rather different from the original one. 

3'1. Two-Zone propagation model. - The modelization of cosmic rays propa- 
gation is a hard task due to the propagation complexity. We use a successfully tested 
model for propagation of cosmic rays [19, 20] which is tuned properly to explain the 
current observations on abundances of nuclear cosmic rays. In general terms, this model 
is divided in two parts: 

• The propagation zone: the region related to the extension of turbulent magnetic 
fields, where cosmic rays propagate in a diffusive regime. Outside this region, we 
suppose cosmic rays propagate in straight lines: this induces cosmic rays leaking 
from the diffusive zone. 

The propagation zone is composed by two cylinders centered at the Galactic cen- 
ter (Figure 3). Both cylinders have a common radius equal to the galactic one 
{R ~ 20 kpc). The thick cylinder has a height of 2L and fills the whole propa- 
gation zone. Its height is constrained by measurements of many species of cosmic 
rays [20]. The second cylinder is a thin disk with height 2hz, where hz = 100 pc. It 
also contains the interstellar medium, cosmic rays sources and interactions among 
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them. 



• The transport equation: it corresponds to a continuity equation for the number 
density of cosmic rays per unit of energy. It contains all physical processes that 
are related to cosmic rays physics, like energy losses, diffusion, re-acceleration and 
cosmic rays sources. 

3'2. Transport Equation for electrons and positrons. - In this case, the 
transport equation is slightly different from the one for nuclei cosmic rays because the 
dominant energy loss process is the interaction with radiation fields (via inverse Compton 
scattering) which is much efficient in these particles than in nuclei cosmic ray. 
Then the transport equation for number density of electron and positrons (tJj) is: 



where D{e) is the diffusion coefficient, b{e) the energy loss term, and q the source term. 
The diffusion term is considered homogeneous in space with a energy dependence: 



where Kq and S are phenomenological parameters inspired by models of the interstellar 
medium based on magneto hydrodynamics. Eq is a normalization energy scale, here fixed 
at the value of 1 GeV. 

Energy losses are related to interactions of electron and positrons with the interstel- 
lar radiation fields, the cosmic microwave background, and the galactic magnetic field. 
Commonly, this term is: 



where te{~ 10^^ s) corresponds to an effective energy loss time calculated via the inverse 
Compton scattering with the radiation fields in the Thomson regime. However, this term 
is no longer valid for energies larger than mg/(e^), where (e^) is the radiation field en- 
ergy, and it has to be corrected considering the Klein-Nishina cross section [21, 22, 23]. 

In Figure 4, we present the energy loss term for the different components of the 
radiation field and those are compared to the standard assumption of the Thomson 
regime (Equation 7). The figure is from our recent work [24] that follows the line of the 
Author's thesis. 

3'3. Transport equation solution. - We solve the differential caution by means 
of Green functions which correspond to the solution for a point-like source: 



(5) 




(6) 




(7) 




(8) 



q{t,x,e) = 6{t - ts)5{e - es)S^{x ~ Xs) , 
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Energy [GeV] 

Figure 4. - Energy loss term 10^''e^^6(e) versus electron (positron) energy. The energy loss for 
different radiation fields models and the standard assumption are shown 



where ts, Ss, and Xg correspond to the injection time, energy, and position. 

The green function is easily obtained through the Fourier transformation method [25, 
1] and it has an analytical form [24]: 

(9) G(t, to, X, X., £, £.) = S{t r,) ^^y;'^*^^ 

5(e) 

^(t tg Tf.) Ggg (X, Xg , £ , e^) , 

where G is the tilded green function: 

~ 1 f [x — xo)^ 

(10) G{x,xo,\d) ^ ^3/2x^ \ 

and it is proportional to the Green function (Ggs) in the time independent case (steady- 
state regime). 

The terms and \d are rich in physical meaning and come naturally from the solution. 
Tc is the cooling time which corresponds to the time lapses for an electron or positron to 
reduce its energy from Ss to e: 

(11) r,(e,e,)= / ' de l/6(e) . 



In the same way, is the diffusion length which gives a scale for the propagation and 
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is defined by: 

(12) A2(£,£,)=4/ deDie)/bie). 



The former two quantities become fundamental to describe the propagation. Also those 
can be calculated for any type of diffusion coefficient and energy loss term allowing us 
to explore new forms, like: deviations from the standard diffusion coefficient and extra 
energy loss processes. 

With the Green function, we are able to find the solution for the transport equation 
for any generic source Q. The solution is found by convoluting the source with the green 
function in the following way: 

(13) tp{t,x,e)= dts dss d^Xs G{t,ts,x,Xs,e,es) X Q{ts,Xs,es) , 



where the integration limits are set to consider only the physical contribution from the 
source Q. 

As we previously said, the model considers that cosmic rays may escape from the 
propagation zone when they reach its boundaries. This means: 



(14) z = ±L\/ ^/x^ +y^' = R =^ '4j{t,x,e) = 0, 

these conditions just affect the tilded Green function (Equation 10). 

There are two standard approaches to impose the boundary conditions: 

• Eigenfunction expansion: The transport equation is solved using a complete set 
of Hclmholtz equation cigenfunctions (xg) that naturally respect the boundary 
conditions. The green function is analytical in some cases and correspond to the 
sum on eigenvalues (17): 

(15) Ghc(x,Xs,Xd) = '^xl{xs)xg{x)cxp (-jg'^Xl 



We use the Fourier-Bessel expansion, where harmonics functions are used for ver- 
tical coordinate and first kind Bessel functions for the radial one. The method is 
very well behaved and fast to calculate when is rather big [26, 27, 1]. 



• Method of Images: Also known as Method of Inversion, it consists on adding counter 
terms or mirror sources that compensate the effect of the original source, preserving 
the boundary condition. 
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For the vertical coordinate, we include an infinite number of mirror sources trans- 
forming the green function into a sum: 

oo 

(16) G'verticalbc(^,^s,Ad) = ^ (-1)" G(f , Xs_„, Ad) , 



n— — oo 



where Xs^n = {^s ,ys ,'2Ln + {—l)"zs) and G is the free space tilded green func- 
tion (Equation 10). This form is very efficient for small values of and complements 
the expansion in harmonic functions. 

The Green function that respects radial boundary conditions is obtained by add 
just one extra counter term: 

(17) Gra,dia.\ hc{x,X s,Xd) = G{x,Xs,Xd) - G(ct, (?s,Ad) 

s Us X . I" y 

where a = (/3a;, fiy, z), a, = (-^, -j,Zs) and = " Ji " ' 

We remark that this Green function works for all values of A^ and does not need 
any kind of special symmetry like the Bessel expansion. 

3'4. Space of Parameters of the model. - The two-zone propagation model 
is also used for studying nuclei cosmic rays. Previous studies on the ratio Boron and 
Carbon (B/C) constraint the space of parameters to a volume which is fully compatible 
with current observations in many other species. We use the allowed space of parameters 
to study positron and electron cosmic rays consistently and systematically with respect 
to cosmic rays propagation (details in [20, 28]). 



4. Positron cosmic rays from dark matter annihilation 

Secondary positrons and electrons are produced in the Milky Way from the interac- 
tion of nuclei cosmic rays on the interstellar gas [29] and are an important tool for the 
comprehension of cosmic-ray propagation [20] . Data on the cosmic positron flux (often 
reported in terms of the positron fraction) have been collected by several experiments 
[30, 31, 32, 33, 34, 35]. 

We point out the HEAT balloon experiment [30] that has mildly indicated a possible 
positron excess at energies larger than 10 GeV with respect to the current - at that 
moment - calculations for the secondary component [29]. In October 2008, the latest 
results of PAMELA experiment [36] have confirmed and extended this feature [37] . 

Different astrophysical contributions to the positron fraction in the 10 GcV region 
have been explored [30] , but only accurate and energy extended data could confirm the 
presence of a bump in the positron fraction and its physical interpretation. Alternatively, 
it has been conjectured that the possible excess of positrons found in the HEAT data 
could be due to annihilation of dark matter in the galactic halo [25, 3]. Although, this 
interpretation is limited by uncertainties in the halo structure and in the cosmic rays 
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Table II. - Dark matter distribution profiles in the Milky Way. 



propagation modeling. 

This section is based on our work [26] . We study the propagation of the positrons 
related to dark matter annihilations in connection with the study of the uncertainties 
due to propagation models compatible with B/C measurements [20]. 



4'1. DARK MATTER ANNIHILATION LIKE SOURCE OF POSITRONS. - According tO the 

various supersymmctric theories, the annihilation of a dark matter particles points to 
the direct production of an electron-positron pair or to the production of many species 
subsequently decaying into photons, neutrinos, hadrons and positrons. 
In our study, we consider four possible annihilation channels which would appear in any 
model of dark matter. The first corresponds to direct production of a e'^e~ pair and it is 
better motivated in theories with universal extra-dimension [41, 42, 43]. We alternatively 
consider annihilations into W^W^ , t^t^ and bb pairs. 



For any annihilation channel, the source term is written as: 

(18) qdmi^i^e) = V {<^v) 



where the 77 is a quantum coefficient which depends on whether the particle is or not its 
own antiparticle. (av) corresponds to the thermally averaged annihilation cross section, 
its value depends on the specific supersymmetric model and it is also constrained by cos- 
mology [44, 45]. We have actually taken here a benchmark value of 2.1 x 10~^^ cm^ sec~^ 
which leads to a relic abundance of ^^^h'^ ^ 0.14 (in agreement with the WMAP obser- 
vations [46, 47]). The dark matter mass value (m^) is still unknown but for the case of 
neutralinos, theoretical arguments as well as the LEP and WMAP results constrain its 
mass to range from a few GeV [48, 49, 50, 51] up to a few TeV. Keeping in mind the 
positron excess, we explore the cases with a dark matter mass of 100 GeV and 500 GcV. 

dn/dse represents the multiplicity distribution of electrons (positrons) per single an- 
nihilation (details in [1]). 

The astronomical ingredient on the source term (Equation 18) is the dark matter 
distribution /9(x) inside the Milky Way halo. We have considered the generic profile: 



(/3-7)/q 



where tq = 8.5 kpc is the galactocentric distance of the Solar System. Note that r 
denotes the radius in spherical coordinates. The local dark matter density has been set 
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equal to pq = 0.3 GeV cm~'^ [52]. We discussed three profiles: an isothermal cored dis- 
tribution [38] for which rg is the radius of the central core, the Navarro, Frenk and White 
profile [39] (hereafter NFW) and Moore's model [40]. However, some studies point out 
to that the central cusp may be less favorable in spiral galaxies and proposed a universal 
dark matter distribution [53, 54]. The NFW and Moore profiles have been numerically 
established thanks to N-body simulations. In the case of the Moore profile, the index 
7 lies between 1 and 1.5 and we have chosen a value of 1.3 (Table II) which is more 
representative. 

The possible presence of dark matter substructures inside these smooth distributions 
enhances the annihilation signal by the so-called boost factor, although the boost factor 
value is still open to debate [55, 56, 57]. It has recently been shown that the boost factor 
due to substructures in the dark matter halo depends on the positron energy and on 
the statistical properties of the dark matter distribution [58]. In addition, it has been 
pointed out that its numerical values is quite modest [59], being of the order of 10-30. 

4'2. THE POSITRON FLUX AND FRACTION. ~ Using the source term (Equation 18), we 
calculate the propagated positrons flux for many cases (Figure 5): 

• Annihilation channels: The nature of dark matter fixes somehow the annihilation 
channel. We observe how channels - like the direct production one - produce 
harder spectra than channels involving quarks because positron are produced with 
energy equal to dark matter particle mass. This is directly related to the multiplic- 
ity distributions. Positrons in quark-antiquark channels come from hadronization 
processes and produce softer spectra with behavior similar to a power law: e~^^'^'^. 
On the other hand, channels involving muons or gauge boson produce harder 
spectra since positrons are produced at earlier stages, in decay chains of the original 
particles, taking a big fraction of the available energy. 

• Propagation uncertainties: Different propagation scenarios have different impact 
on positron generated by dark matter annihilation. In Figure 5, we present the 
uncertainty band associated to the B/C analysis. We notice that the size of the 
band depends on the annihilation channel. The direct production case presents 
the largest uncertainty band at low energies because the low energy flux is the 
result of far-away propagated positrons and not from the very local ones. This 
also explains why at energies closer to 100 GeV, the uncertainty band is smaller 
converges to an unique curve. Other channels arc less affected because very local 
produced positrons are not affected by the propagation parameters. 

Let us stress that the positron flux obtained from annihilation of dark matter is not 
the only one. We need to consider other astrophysical components. Due to the nature of 
astrophysical processes, positrons are dominated by a secondary component, i.e. those 
are created from the interaction of nuclei cosmic rays with the interstellar gas. 

To study the behavior of the positron signal, we include the secondary positron com- 
ponent and the electron flux from parametrized fluxes [25]. 

In Figure 6 we present the effect of the annihilation channel and propagation uncertainties 
on the positron fraction. We observe how channels with harder spectra arc more suitable 
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Figure 5. - Positron flux versus the positron energy e, for a dark matter particle with a 

mass of 100 GeV and for a NFW profile (Table II). The four panels refer to different annihilation 
final states: direct e^e" production (top left), 6b (top right), W^W' (bottom left) and t^t~ 
(bottom right). In each panel, the thick solid [red] curve refers to the best-fit choice (MED) 
of the astrophysical parameters. The upper [blue] and lower [green] thin solid lines correspond 
respectively to the astrophysical configurations which provide here the maximal (Ml) and min- 
imal (M2) fiux. The numerical values of these configuration are defined in [26]. The colored 
[yellow] area features the total uncertainty band arising from positron propagation. 



to explain the positron excess. The hh case is less favorable because most of the positrons 

are at low energy, making impossible to reproduce the observations. 

The uncertainties of propagation are sizable with respect to the signal from dark matter, 

however the impact is not enough to destroy some features arising from the annihilation 

channels. 



5. Secondary positron flux at Earth 

The secondary component is the result of the interaction of cosmic rays nuclei with 
the interstellar gas composed mainly by hydrogen and helium. We model the interstellar 
gas assuming a homogeneous disk with radius equal to the galactic one and thickness of 
200 pc (details in [1]). 
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Figure 6. - Positron fraction /{e~ + e"*") versus the positron detection energy e. Notations 
are as in Figure 5. In each panel, the thin [brown] solid line stands for the background [25, 
29] whereas the thick solid [red] curve refers to the total positron fiux where the signal is 
calculated with the best-fit choice (MED) of the astrophysical parameters. Experimental data 
from HEAT [30], AMS [32, 33], CAPRICE [34] and MASS [35] are also plotted. 



The source term of secondary electrons and positrons coming from interaction of 
proton cosmic rays with hydrogen is: 



(20) 



(7ph(x, Ee) = 47r 71h(x) j dsp 3>p(x, Sp) 



dee 



(Spj £e) J 



where nn is the hydrogen number density, $p is the proton flux (details in [1]), and 
(rja) represents the inclusive cross section of the process p + p — > e^(e~) + X discussed 
and calculated in Section 2. Let us remained that the inclusive cross section, in our case, 
comes from the invariant cross section parameterizations of Badhwar et al. [9], Tan and 
Ng [12], and Kamae et al. [14]. 



Let us stress that the complete positron and electron source term is the sum of all 
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Figure 7. - Interstellar electron and positron fluxes e'^^e for the MED set versus energy. The 
curves correspond to fluxes calculated from Nuclei CR interactions with the ISM. Each curve 
represents a different nuclear cross section and Nuclei CR flux parameterization. Also the 
Strong et al. [25] flux parameterizations are shown. The uncertainty band related to those 
parameterization is plotted (yellow band) as well. 



contributions from proton, alpha particles, hydrogen and helium: 

(21) (?fun(x,£e)=47r ^ ^ „^(x) / de, $,(£,) ^^^(e„ £e) , 

2— H,Hc j—p^a 

where the inclusive cross sections for processes p+He, a+H and a+He are estimated by 
scaling proton-hydrogen cross section: 

where Aj is the mass number of the incident particle and s/ are scaling factors [60, 61]. 



5'1. SECONDARY POSITRON AND POSITRON FRACTION. - The propagation of sec- 
ondary positrons is realized according to the two-zone propagation model (section Section 3). 
Similar to the case of annihilation of dark matter, we study the effects of uncertainties 
related to nuclear physics and propagation. In Figure 7, we present secondary electrons 
and positrons obtained from the three nuclear parameterization: Badhwar et al. [9] , Tan 
and Ng [12] and Kamae et al. [14]. We observe that at high energy (> lOGeV) all cases 
converge because at that energies all three parameterizations converge due to most of 
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Figure 8. - Secondary positron flux £^'^$e+ versus positron energy. Positron fluxes were calcu- 
lated using proton and alpha CR fluxes from Donato et al. [62] and Shikaze et al. [8] with the 
Tan et al. cross section parameterization [12]. 



experimental data is in this range. At low energy, we observe how differences on the nu- 
clear parameterization affect the propagated positrons and electrons. Nevertheless, this 
effect is not so important because it occurs in the energy range where solar modulation 
dominates. 

Moreover, we can not forget the presence of the uncertainties related to propagation. 
Figure 8 shows the secondary positron flux versus positron energy. We consider the whole 
space of parameter compatible with the B/C analysis [20] , obtaining that the estimated 
flux encompasses the data. This confirms the compatibility among electron-positron 
cosmic ray propagation with nuclei cosmic rays case from the propagation point of view. 




Figure 9. - Positron fraction as a function of the positron energy, for a soft and hard elec- 
tron spectrum. Data are taken from CAPRICE [34], HEAT [30], AMS [33], MASS [35] and 
PAMELA [37]. 
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The natural extension of the analysis is to calculate the positron fraction. Using 
the calculated secondary flux and the total electron flux estimated from observation like 
AMS, we observe how the positron fraction is highly sensitive to small variation in the 
electron flux (Figure 9). We find that harder spectra makes stronger the evidence of 
positron excess, however a softer electron spectra makes possible to explain the raise in 
the faction using mainly secondary positrons. 

This study points towards a more detailed study on the electron flux which has been 
recently addressed in Rcf. [24]. 



6. Conclusions 

During the last years, many of the latest cosmic rays experiments have shown very 
interesting results that are pushing to new frontiers the knowledge about the galactic 
environment and the origin of the GeV-TeV cosmic rays. And of course, the case of 
electron and positron cosmic rays is not an exception. 

Dark matter annihilation is a very exciting possibility to explain the positron excess, 
which is not sufficiently explained by contributions from secondary positrons. In our 
works, we studied propagation uncertainties associated to the analysis made on cosmic 
rays nuclei (B/C). These uncertainties affect considerably the signal coming from anni- 
hilation of dark matter, however, most of the characteristics related to the annihilation 
were not significantly modified. This fact is promising for further analysis on the signal. 

To go deeper in the analysis, we studied the secondary production of positrons. We 
considered uncertainties related to nuclear physics in addition to the ones related to prop- 
agation. One of the first results was that secondary positron production is compatible 
with current measurements and with the B/C analysis, which is remarkable considering 
that we proposed that cosmic ray propagation is common for all the species. 

In both studies, the positron fraction was analyzed. Dark matter annihilation sce- 
nario is able to explain it, although, the lack of precision from the theoretical point of 
view makes hard to identify an exotic component present in it. Moreover, the positron 
fraction is very sensitive to variations in the electron fiux. 

We stress the necessity to re-estimate secondary and primary electron component, 
and to consider already known sources, especially pulsars, in order to discard/confirm 
possible presence of an undiscovered component, like the case of dark matter annihila- 
tions. 
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